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ABSTRACT 


yntafriy  beaea th  the  suifucu  uf  g'fhild  will  uure’i 
urnhs  4m  Wfcwwd  UU  BWTtMeiy*' The  detection  of  t£s  wake  le  usually  very  difficult 
because  of  its  small  amplitude,  and  because  of  the  natural  waves  which  are  usually 
present.  It  is  shown  that  the  three  dimensional  spectra  of  the  wake  and  of  the  natural 
waves  are  of  different  forms,  and  that  this  difference  might  be  used  to  design  a  filter 
to  enhance  the  wake  and  attenuate  the  natural  waves.  Some  of  the  practical  limitations 
on  such  a  filter  are  also  discussed. 


f 


PIBMRI- 1037-62 


TABLE  OF  CONTENTS 

Page 

Introduction  1 

Conditions  of  the  problem  1 

Spectrum  of  gravity  wave  noise  2 

Signal  spectrum:  wake  of  a  moving  submerged  object  6 

Filtering  in  two  dimensions  7 

References  13 

Appendix  1  -  Wake  of  Submerged  Cylinder  A- 1 

Appendix  2  -  Width  of  the  Spectrum  of  Gravity  Waves  A-2.  1 

Finite  depth  A-2.  1 

Capillarity  A-2.  1 

Finite  Amplitude  A-2.  1 

Nonlinear  superposition 


A-2.  2 


PIBMRI- 1037-62 


FILTERS  FOR  DETECTING  MOVING  UNDERWATER  OBJECTS 


by 


Arthur  E.  Laemmel 


L  Introduction  The  problem  of  detecting  the  surface  wake  of  an  object  moving 
beneath  a  fluid  surface  can  be  considered  as  analogous  to  the  familiar  problem  of 
detecting  an  electrical  signal.  In  both  cases  the  desired  effect  is  partially  hidden  by 
a  random  disturbance:  surface  gravity  waves  in  the  case  of  the  wake,  and  antenna  or 
tube  noise  in  the  case  of  the  electrical  signal.  Several  types  of  filter  for  separating 
the  desired  effect  or  signal  from  the  random  disturbance  or  noise  have  been  extensively 
studied  in  the  electrical  case,  for  example  the  filters  of  North^  and  Dwork,  and 
Wiener^  Some  extensions  have  been  made  to  more  than  one  dimension  in  the  case 
of  optics,  infra  red,  and  radar;  but  these  have  generally  been  applied  to  several  spatial 
dimensions  without  considering  the  time  dimension  simultaneously.  The  filters  ana¬ 
lyzed  below  are  similar  to  the  MTI  (moving  target  indicator)  radars,  but  the  special 
properties  of  water  waves  introduce  several  novel  and  promising  opportunities  for 
improvement. 

2.  Conditions  of  the, problem  Let  a  rectangular  coordinate  system  be  established 
with  the  origin  at  the  undisturbed  surface  of  the  fluid  with  z  going  upward  and  x  and  y 
going  horizontally.  Let  A(x,  y,  t)  be  the  wake  of  the  moving  object  and  B(x,  y,  t)  be  the 
random  disturbance  s  on  the  surface  of  the  fluid.  The  variables  A  and  B  can  represent 
the  value  of  z  at  the  surface,  radar  reflectivity  etc.  ;  the  only  restriction  being  that 
the  noise  B  is  linearly  added  to  the  signal  A  before  being  received  by  the  detector. 

A  three  dimensional  spectrum  will  be  defined  for  A  and  called  a,  and  according  to  the 
usual  Fourier  transform  theory,  A  can  then  in  turn  be  represented  in  terms  of  a: 

ao 

a(u#w  ,£)=  J  J  J  A(x,  y,t)  e~2  *  ^xu  +  ^  +  dxdydt  (1) 

-oo 

oo 

A{x, y, t)  =  /  /  /  a(u,w,f).  2'i<x“  t’rw+,£> 

-00 


du  dw  df 


(2) 
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Analogous-  expressions  can  be  written  representing  B  in  terms  of  its  spectrum  p.  Note 
that  f  is  the  frequency  familiar  to  electrical  engineers  and  is  measured  in  cycles  per 
second,  while  u  and  w  are  the  spatial  frequencies  which  are  beginning  to  be  common 
in  optical  engineering  and  are  measured  in  cycles  per  foot  or  other  unit  of  length. 

Let  y  (u,  w,  f )  be  the  frequency  function  of  a  linear  filter  which  acts  on  the  signal  as 
follows: 


a(u,  w,  f)  =  y  (u,  w,  f)  a(u,  w,  f) 


(3) 


The  representation  of  the  filtering  in  the  space-time  domain  can  be  shown  to  be  the 
following  convolution: 

00 

X(|,4,t  )=//  /  r  (ft -x,  4-y,  T  -  t)  A(x,y,t)  dx'dy  dt  (4) 

-oo 


where  r  is  defined  analogously  to  A  in  Eq  2.  The  function  r  can  be  regarded  as  the 
filter  response  to  a  unit  impulse  at  the  origin  x  =  0,  y  =  0  and  at  time  t  *  0.  Note 
that  a  bar  is  used  to  denote  that  the  corresponding  function  has  been  filtered,  and  15 
and  "{Twill  be  similarly  defined.  The  object  of  the  filter  can  be  very  generally  stated 
as  this:  to  reduce  t^e  function  15  to  as  small  a  value  as  possible  while  keeping  A  large 
at  somewhere  near  the  target' s  position.  The  filter  of  North  and  Dwork  maximises 
the  peak  value  of  A  without  regard  to  whether  or  not  A  corresponds  to  A  in  shape. 
Wiener' s  filter  attempts  to  make  A  +  B  as  much  like  A  in  shape  as  possible.  The 
former  specifies  the  phase  and  amplitude  of  the  filter  spectrum,  the  later  only  the 
power  spectrum.  These  and  several  other  types  of  filter  have  one  simple  feature  in 
common,  i.  e.  they  pass  regions  of  the  spectral  plane  where  the  signal  is  large  and 
attenuated  regions  where  the  noise  is  large.  This  effect  is  intuitively  obvious;  the 
more  formal  theories  are  only  necessary,  where  the  spectra  of  signal  and  noise  over- 
lap,  to  give  the  relative  weighting  in  such  overlap  regions.  The  greater  the  overlap 
(i.  e.  the  similarity  in  shape  of  signal  and  noise  spectra)  the  more  one  must  rely  on 
the  formulas,  but  also  the  less  likely  that  a  successful  separation  of  signal  and  noise 
can  take  place. 

3.  Spectrum  of  gravity  wave  noise  Some  simple  examples  of  noise  and  signal  will  be 
examined  to  give  a  feeling  for  the  relative  disjointness  of  the  spectra.  The  physical 
processes  are  most  completely  understood  for  the  case  where  the  variables  A  and  B 
represent  the  vertical  displacement  of  the  surface  of  the  fluid.  If  these  displacements 
are  small,  the  wake  of  a  moving  submerged  object  will  be  linearly  super  imposed  on 
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the  wavea  normally  present  on  the  surface  of  the  fluid.  An  elementary  wave  on  the 

(4  5 ) 

surface  of  a  fluid  of  constant  depth  h  is  '  ’  ' 

b(x,y,  t)  =  cos  (  u>t  -  kx  cos  0  -  ky  sin  0  +  i|i)  (5) 

where  «  =  +>/g  h  tanh  k  h  (fc) 

i|/=  a  phase  constant 
0  =  ‘angle  of  propagation  from  x  axis: 
g  a  acceleration  of  the  Earth' s  gravity 


2 

The  wavelength  is  seen  to  be  -£■  (=  X).  The  frequency  of  oscillation  at  a  fixed  point 
is  cycles  per  second.  The  phase  velocity  of  the  wave  is 

C  IT 

c  =  TT=  Unh 


It  will  be  noted  that  although  to  and  k  can  assume  any  positive  value  between  0  and  go  , 
the  value  of  c  ranges  between  0  and  v§h  ,  If  the  wavelength  is  long  enough  for  c  to 
be  near  its  upper  bound  the  wave  is  called  a  tidal  wave.  In  the  present  problem  the 
opposite  case  is  of  more  interest:  If  the  wavelength  is  small  compared  to  the  depth  the 
velocity  is  approximately^/^- —  , 
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The  most  seneral  function  describing  the  free  vibrations  of  the  fluid  surface  is 
obtained  by  combining  many  components  of  the  type  represented  by  Eq.  5. 


B{x,  y,  t)«y  p  cos  (w  t  -  k  x  cos  0  -  k  y  sin  0V  +  >|)v. ) 

“V  v 


(8) 


where 


gk„  tanh  h  k^ 


(9) 


A  more  general  representation  is  obtained  if  this  summation  is  replaced  by  an  integral 
and  if  the  substitutions 


2  ir  u  -  -kv  cob  0j, 
2  ir  w  =  -k^  sin  0„ 


(10) 


are  made: 


oo  oo 


(ID 


B(x,y,t)«  J  J  p(u,  w)  cos  j[to  (u,  w)t  +  2  w  ux  +  2  w  wy  +  t|i(u,  w)]  iudw 


-oo  *oo 


where  w(u,  w)  *^2*  g^u2  +  w2)  tanh  2  w  h^  +  w  ) 


(12) 


This  can  be  rewritten  in  the  following  exponential  forms: 

2  ir  i  (ux  +  wy  +  ft) 


oo  oo  oo 


B(x,  y,  t)  *  /  /  /  P(u,w,f) 


du  dw  df  (13) 


-oo  -oo  -oo 


-i»|<(-u,  -w) 


(14) 


r, .  w  (u»  w)  I 

■L** — rs  j 
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oo 

and  6(f)*  0  (f  /  0),  J  6(f)  df  =  1  (impulse  functions) 
•oo 


and  where  to  (  )  is  defined  in  Eq.  12. 

Sereral  very  interesting  properties  of  the  spectrum  of  the  wave  "noise"  'are 
at  once  apparent: 


1)  p  (-u,  -w,  -f)  *  p*  (u,w,f) 

e 

2)  p  (u,  w,  f )  =  0  except  on  the  surface  f  *  • 

(see  Eq.  12) 

3)  If  all  of  the  original  wave  components  are  travelling  in  the 
in  the  same  direction,  say  at  an  angle  6o  with  the  x  axis,  then  p(u,  w,  f)  *  0  except 
along  the  line  defined  by  the  equations: 


U  *  or  COB  e 


w  =  it  sin  0 


f  » 


■yjj-  w  (tr  cos  0o,  v  sin  0Q)  0<  v 
w  (v  cos  0Q  ,  tr  sin  0Q)  0  <  0 


see  Eq  12 
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4,  Signal  spectrum:  wake  of  a  moving  submerged  object  The  simplest  case  to  con¬ 
sider  will  be  that  in  which  the  object  sets  up  a  surface  distrubance  which  travels  with 
it  without  changing  shape.  If  the  object  is  travelling  with  a  velocity  c  in  a  direction 
making  an  angle  0  with  the  positive  x  axis,  then  the  signal  function  can  be  written 


A(5,4  -  ct) 


where 


x  sin  0  -  y  cos  0+1 
x  cos  0  +  y  sin  0 


The  spectrum  can  now  be  calculated  by  Eq.  (1): 


oo  oo  oo 
u(u,  w,  f )  *  f  f  f 
-oo  -oo  -oo 


-2  v  i(ux  +  wy  +  ft) 

A(e,  C-ct)e 


dx  dy  dt 


(15) 


a(ut  w,  f)  *  (X  (u  sin  0  -  w  cos  0  ,  u  cos  0  +  w  sin  0) 


(16) 


6(f  +  u c  cos  0  +  w  c  sin  0)e 
oo  oo 


-2  ir  il(w  cos  0  -  u  sin  0) 


n  2  w  i(£  u  +  £  w) 

where  U  (u,  w  )=  /  /  A(| , ;  )e  ds  d; 


-GO  -00 
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The  following  interesting  properites  may  now  be  deduced: 

$ 

1)  q(-u, -w,  -f)  sa  (u,w,  f) 

2)  a(u,  w,  f )  =  0,  except  on  the  (plane)  surface  f  =  -uc  cos  6  -  wc  sin  0 

3)  The  intercept  of  the  above  plane  with  the  (u,  w)  plane  is  the  line 
u  =  -*tan  0 


Note  particularly  that  this  line  does  not  depend  on  the  distance  of  the  object' s  path 
from  the  origin  (I).  The  angle  between  the  two  planes  depends  only  on  the  magnitude  of 
the  velocity  and  is  fact  tan-  *  c  . 

Thus,  the  first  order  theory  shows  that  the  signal  spectrum  exists  only  on  plane 
in  the  spectral  domain  and  the  noise  spectrum  exists  only  on  a  paraboloid  of  revolution 
(for  infinite  depth).  Offhand,  this  sounds  almost  ideal  from  the  viewpoint  of  filtering. 
Unfortunately,  the  plane  and  the  paraboloid  intersect  and  most  of  the  signal  is  near  this 
curve  of  intersection. 

5.  Filtering  in  two  dimensions  Most  of  the  principles  underlying  the  proposed  methods 
of  filtering  can  be  illustrated  in  two  dimensions,  thus  avoiding  temporarily  some  of  the 
mathematical  complexity  of  the  three  dimensional  case.  It  is  not  suggested  that  the  two 
dimensional  case  indicates  anything  about  the  feasibility  of  detection  in  three  dimensions. 

Assume,  as  before,  a  coordinate  system  with  its  origin  at  the  undisturbed  level 
of  the  fluid,  z  going  upward,  and  with  all  boundaries  and  variables  independent  of  y. 

If  the  depth  is  infinite,  the  "noise”  will  be  free  gravity  waves  of  the  form  (see  Eq.  8) 
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The  convention  will  be  adapted  that  the  Wavelength  will  always  be  positive,  but  that 
can  be  positive  (waves  going  in  the  -  x  direction)  or  can  be  negative  (waves  going 
in  the  1-  x  direction).  The  spectrum  of  a  general  function  of  x  and  t  can  be  represented 
quite  generally  by  drawing  "  contour”  lines  of  constant  amplitude  and  constant  phase 
in  the  plane  of  one  spatial  frequency  (cycles  per  foot)  and  are  temporal  frequency 
(cycles  per  second).  This  spectral  plane  should  not  be  confused  with  the  more  usual 
complex  frequency  plane  where  the  temporal  frequancy  is  given  an  imaginary  part 
(makes  per  second).  In  the  present  case  the  spectrum  vanishes  almost  everywhere  in 
the  spectral  plane  except  along  the  lines  specified  by  Eq.  18.  Along  these  lines  Dirac 
delta  functions  are  distributed,  as  in  Eq.  14. 
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The  signal  to  be  detected  in  the  present  example  will  be  taken  as  the  wake  o£ 
a  circular  cylinder  of  radius  b  moving  with  a  constant  velocity  c  in  the  -  x  direction 
at  a  constant  depth  s  =  -  D.  A  modification  of  the  method  used  by  Lamb  is  given  in 
Appendix  1  and  results  in: 


A(x,  t) 


,.2  2 
2b  c 


g 


x  +  ct)2  -  D2 
(x  +  ct)Z  +DZ]  2 


go 


it  (x  +  c  t)  sii 


(19) 


where  r 


<4)  = 


{ 


o 

l 


4  <  o 
o<  4 


The  spectrum  can  be  obtained  from  an  intermediate  result  in  Lamb*  s  method,  since 
he  actually  satisfied  the  surface  boundary  condition  in  terms  of  the  spectral  function. 
From  Eq.  <  12  in  appendix  1: 


oo  -kd 

A(x,  t)=  2  b 2  J  -^-5 - - -  cos  (kx  +  kct)  dk 


i  x-4- 


(20) 


00  00 


A(x,t)=  J  J  bZ  6(w -kc) 


Ikl  e”  k  D  i(hx+w  t) 

L^-L- = -  e  du  dk  (21) 


-00  -00 


lkl  — $— 


Thus 


a(u, f)  =  4x2  b2  6(f-uc)  iHil. 


-2  w  u  D 


u  - 


2  IT  c 


(22) 


This  vanishes  everywhere  in  the  (u,  f)  plane  except  along  the  line  f  =  u  c. 
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The  filter  which  is  to  separate  signal  from  noise  can  be  designed  in  several 
ways:  reject  frequencies  which  represent  free  surface  waves,pass  only  those  fre¬ 
quencies  representing  a  moving  object,  or  to  pass  or  attenuate  any  particular  region 
according  to  its  effect  on  the  signal  to  noise  ratio  or  probability  of  detection.  The  first 
is  perhaps  simplest  and  does  not  require  knowledge  of  the  object' s  velocity.  To  see 
the  effect  of  such  a  filter  on  the  desired  wake,  Eq.  21  must  be  evaluated  over  the  whole 
plane  except  for  a  band 


-  A  <  <*>  <  t  f  A  (23) 

where  the  -  is  to  be  interpreted  as  the  jgame  sign  ask. 

This  is  most  easily  done  by  subtracting  the  contribution  from  this  band: 


oo  oo 

Aq(x  ,t)  -  J  J  yQ  (k,w)bZ  6(io-kc) 
-oo  -oo 


|k|  e~  ^  i(kx  +  u>t) 

w  --V 

c 


(24) 


where 


Y  Q  (k,  u>) 


1  Tgk  -  A  c  uX  Jifik  +  a 
1-  m  -  A<  <*>  <  lkl  +  A 
>  otherwise 


0  <  k 
k  <  0 


oo 

AQ(x,t)  =  /  Y0(k,kc)bZ  _liile 


-|k|D  ik(x  +  ct) 


-oo 


W  --V 


dk 


(25) 


This  can  be  approximated  as  follows  for  small  a 


,  A1  .  -kD 

A  (x,t)  «  2b  f  — - - -  cos  k(x +  ct)dk  +  2b6  f 

O  J  K  J 


"V  + 


o  c  i  ^  -kD 

2  r  ke 


«■-+ 


■“-V 


cos  k(xi-ct)dk 


The  bandwidths  a  j  and  ^  are  functions  of  a  and  other  parameters.  They  approach 
zero.  Their  calculation  is  discussed  below. 
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Neglecting  terms  involving  powers  of  and  ^  higher  than  the  first: 


Aq  (x,  t){ 


4b 


e  c  Si  ^Ix  +  ctlj  8in  ■+  |x  +  ct| 


(26) 


The  value  of  the  integral  over  all  space  is 


Aco  <*•*> 


-  2bZc2  (x  +  ct)2  -  D2 


6 


x+ct)2  +  D2] 


2  "-*? 

— e  C  sin |  x -ret  | 
c‘  c 


(27) 


The  wake  function  after  filtering  by  the  band  elimination  filter  is  therefore  the  differ¬ 
ence  between  these  two: 


X(x,t)=i^£-2  f«.+  *4-P-L 
g  [(x+ctr+D*] 


gD 

C  {^-F  Si  U|x  +  ct|]}  sin&y 


x  +  ct 


(28) 


The  disturbance  represented  by  Eq.  28  is  seen  to  be  composed  of  two  parts: 

2b2c2 

1)  a  symmetrical  depression  of  amplitude - «—  ft.  centered  over  the  cylinder  and 

gD 


of  width  2  D  ft.  or 


2  D 


sec,  and  2)  a  packet  of  sine  waves  of  amplitude 


■f 


-»D 
c 


'¥■  2 

“  c  c  2 

,  of  period  2  w  — — —  ft.  or  2  ir  —  sec. ,  and  of  duration -  ft. 

*  8  AZ 


2 

or  — —  sec.  The  later  component  is  also  symmetrically  centered  over  the  cylinder, 

^2C 

but  for  certain  types  of  filter  it  may  be  behind  the  cylinder  because  of  physical  reali¬ 
zability  conditions  on  the  filter.  For  various  reasons,  some  of  which  are  analyzed  in 
Appendix  2,  the  noise  spectrum  is  not  confined  to  a  line  w  =  •/$£  .  If  the  noise 
(surface  gravity  waves)  is  to  be  eliminates,  then  a  band  of  frequencies,  or  several 
bands  of  frequencies,  must  be  rejected  by  the  filter.  One  of  the  principal  reasons 
for  the  zero  spectral  width  is  the  deviation  from  simple  linear  waves  caused  by 
finite  wave  amplitude.  Thus,  the  higher  the  waves,  the  wider  the  spectral  rejection 
band,  and  the  shorter  the  wave  packet  in  Eq.  28  will  be.  A  limit  is  reached  when  the 
wave  packet  is  reduced  to  a  single  loop. 
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Some  idea  of  the  size  of  the  various  parameters  can  be  obtained  from  the 
following  table: 


c 

knots 

X 

feet 

t 

sec. 

2a 

feet 

2a2 

feet 

2/ 

a2 

feet 

cycles 

2 

2.23 

.  66 

0.  2 

0.  028 

17.  8 

8 

5 

13.  9 

1.  65 

1.  25 

0.  18 

111 

8 

10 

55.  7 

3.  30 

5 

0.  71 

446 

8 

20 

223 

6.  60 

10 

0.  71 

7140 

32 

50 

1390 

16.  5 

10 

0.  11 

1,  730,  000 

1240 

The  second  and  third  columns  are  the  wavelength  and  period  respectively  of  the 
oscillitory  part  of  the  wake  of  an  object  moving  at  the  indicated  speed.  These  two 
columns  also  characterize  the  gravity  waves  which  will  cause  the  most  trouble  in 
detection  of  the  moving  object.  The  fourth  column  shows  an  assumed  height  for  the 
gravity  waves  naturally  present.  The  height  is  double  the  amplitude  of  fundamental 
component,  but  this  is  less  than  the  peak  height  usually  measured.  The  values  shown 
in  the  fourth  column  were  arrived  at  as  follows:  The  maximum  peak  height  of  a  wave 
is  known  from  hydrodynamic  theory  to  be  about  1/7  of  its  length^’  For  shorter 
wavelengths,  the  double  amplitude  of  the  fundamental  component  is  taken  as  1/  5.  57 
of  its  length.  For  larger  wavelengths,  common  experience  shows  that  the  amplitude 
does  not  increase  indefinitely  with  wavelength.  Hence,  for  X  =  55.  7  ft.  and  above 
fundamental  double  amplitude  is  kept  at  10  ft.  The  fifth  column  shows  the  double 
amplitude  of  the  second  harmonic  of  the  gravity  waves,  computed  according  to  Eq.  A3 
of  Appendix  2.  The  sixth  column  shows  the  number  of  cycles  in  the  filtered  oscillitory 
wake  when  the  gravity  wave  rejection  filter  is  adjusted  to  have  a  bandwidth  given  by 
Eq.  A 14. 
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Appendix  1  -  Wake  of  Submerged  Cylinder 


(4) 

Refer  to  Lamb  Section  247,  p.  410,  Everything  checks  up  to  Eq.  12,  and  this  is 
correct  formally,  but  it  does  not  remove  the  dominant  contribution  of  the  integral. 
Retain  the  first  form  of  Eq.  12  ^2  ™  e~^  cos  kx 


/ 


k  -  k 


dk 


(L12V 


and  transform  the  variable  in  this  just  as  Lamb  does.  Instead  of  his  Eq.  14,  this 
gives: 


1=  -  2 


,2  -  xf 

it  k  b  e 


sin  k  x  +  2  b 


/ 


K  COS 


mf  - 
~T 


+  m 


m  sin  m  f 
2 - 


-mx  . 
me  dm 


Expanding  the  denominator  in  a  power  series: 


~t~ — T  "  “ 

k  +  m  x 


1  2.  m 

- m  + 


Integrating  term  by  term,  arranged  according  to  ascending  powers  of  m: 


7  -xf  2  2  ,2 

t  =  -  2  *■  «  b  e  sin  k  x  + -  — - 

(x2  4  f2)2 


This  is  for  x  >  0 ,  for  x  <  0  note  that  £  is  an  even  function  of  x  in  Eq.  12.  Add  a  free 
gravity  wave,  just  as  Lamb  does,  to  satisfy  the  upstream  boundary  condition.  Instead 
of  his  Eq.  17,  there  results: 


&(x) 


2b2c2 


-f2 


(T7 


TP 


-  e  f 


s  in 


0<x 


V  (L17 


;(x)  « 


-.2  2 

2b  c 


g 


xi-f2 

(xZ+fZ)Z 


x  <  0 
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Appendix  2  -  Width  of  the  Spectrum  of  Gravity  Waves 

If  the  spectrum  of  the  gravity  wave  "noise"  is  really  confined  to  a  line  (one 
space,  one  time  dimension)  or  to  a  surface  (2  space  one  time  dimensions)  in  the 
spectral  space,  then  the  "noise"  can  be  reduced  indefinitely  allowing  extremely  small 
wakes  to  be  detected.  As  was  pointed  out  in  Section  5  above,  the  sensitivity  of  a 
detector  using  the  filtering  method  proposed  in  this  report  depends  directly  on  how 
narrow  the  rejection  band  can  be  made.  Therefor,  some  estimates  will  be  made  here 
on  how  wide  a  spectral  region  the  actual  waves  occuppy.  Note  that  in  most  cases  the 
spectrum  is  still  a  line  (or  surface)  but  shifted  from  the  normal  position.  This  allows 
the  possibility  that  a  filter  could  be  adopted  to  the  new  situation,  provided  that  the 
new  conditions  apply  consistantly.  For  example,  a  wave  having  an  amplitude  equal 
to  1/  5  its  wavelength  would  merely  have  its  infinite smal-amplitude  spectrum  displaced*, 
but  it  is  then  unreasonable  to  assume  that  smaller  waves  would  not  also  be  present 
similtaneously. 

Finite  depth  From  Eq.  6,  the  frequency  of  a  single  component  plane  gravity 
wave  can  be  seen  to  be 


-  i/gk  tanh  h  k  «v^k  (1-e 


(A1 


Capillarity  From  Lamb,  p.  459  Eq.  3,  a  rearrangement  gives 


<*>  =Vsk  + 


"k"^  T  k^  T 

- L.  «y^k  (i  + - -i  =v£h  (n-  l-  )  (A2 

p  2  g  p  653 X 


where  T,  =  surface  tension  =  74  ^X5.?s. .  =  0.  163 
1  cm 


lb. 

sec 


p  =  density  =  64  ^Vft^ 


Finite  amplitude  Refer  tu  Lamb,  p  417.  The  vertical  displacement  of  the 
surface,  disregarding  the  "  DC"  level,  is  given  by  his  Eq.  3 


,2  ,.2  3 

£  (x,  t)  =  a  cos  k(x  +  ct)  + -  cos  2k  (x  +  ct)  + - - —  cos  3k  (x  +  ct) 

2  8 


(A3 


Harmonics  could  also  be  added,  these  also  along  lines  or  surfaces. 
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Note  that  the  wave  is  considered  to  be  movine  here  instead  of  the  whole  body  of  water. 
His  Eq,  6  gives 


=V^k  U  k2  aZ)«ygk 


(A4 


A  *  g* 


.2  2 

k  a 


In  this  case  another  effect  is  also  observed,  namely  that  not  only  is  the  fundamental 
frequency  shifted  slightly  upward,  but  also  distinct  harmonic  terms  are  added.  The 
time  and  space  frequencies  of  these  harmonics  can  be  expressed  in  terms  of  the 
fundamental  k  as  follows: 


Also,  the  harmonic  time  frequencies  can  be  expressed  as  functions  of  the  harmonic 
space  frequencies: 


Nonlinear  superposition  The  simple  theory  of  gravity  wave  noise  breaks  down 
unless  the  waves  are  of  infinitesimal  amplitude.  One  such  breakdown  is  considered 
above  under  N  finite  amplitude,  n  and  this  is  analogous  to  harmonic  distortion  in 
electrical  systems.  Another  breakdown  which  occurs  in  analogous  to  intermodulation 
distortion.  Under  finite  amplitude  above  only  waves  which  are  periodic  were  consider¬ 
ed.  If  two  waves  of  different  periods  are  super inposed  they  act  independently  of  each 
other  if  they  are  both  of  infinitesimal  amplitude,  i.  e.  the  displacement  of  the  surface 
with  both  waves  present  is  the  sum  of  the  displacements  when  the  two  waves  exist 
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separately.  The  spectra  also  add  linearly,  and  therefor  the  gravity  wave  noise 
spectrum  is  confined  to  the  parabolic  line  mentioned  above. 


Consider  now  the  case  where  two  waves  are  present  at  the  same  time,  one  of 
much  longer  wavelength  than  the  other.  If  x  is  the  distance  along  the  surface  and  £ 
the  elevation  of  the  surface  the  component  with  the  longer  wavelength  can  be  represent¬ 
ed  as 


t  j  (x, t) 


(A7 


Now  if  another  component  is  present  with  a  wavelength  («  X  j),  then  it  is  subject  to 
downward  body  forces  from  not  only  gravity,  but  also  from  the  vertical  acceleration  of 
the  water' s  surface  due  to  the  first  component.  The  effective  acceleration 'of  gravity 
for  the  second  component  is  thus: 

a2 3 

g  (*,  t)  =  g  +  -2-y  C  I  4)  (A8 

at* 


T (x,  t)  s  g  -  2  ir  g  cos^q-  x  +J-2-^q—  t  (A9 

It  is  assumed  that  the  second  component  will  have  little  effect  on  the  first.  The 
problem  of  finding  a  free  surface  wave  under  the  influence  of  a  sinusoidally  varying 
gravity  force  is  analogous  to  that  of  obtaining  frequency  modulation  by  varying  the 
capacitance  in  an  oscillator  circuit..  The  solution  cannot  be  obtained  by  merely  sub  - 
stituting  a  variable  g  function  in  Eq.  A8. 

The  method  of  solution  for  the  smaller  wave  will  be  a  simple  modification  of 
that  given  by  Lamb^  and  Stoker  First,  solve 

2  2 

=  0  (A10 

2  x  2  z 


in  the  interior  of  the  liquid  subject  to  the  following  boundary  conditions: 


=  0  for  z  =  -  oo 

3  z 


±±.  +  -g94_=0  for  z  =  0 


(All 

(A  12 
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Finally)  the  surface  displacement  caused  by  the  second  component  is  given  in  terms 
of  4  by:' 


w-Mf 


E  =  0 


(A  13 


Assume  a  solution  of  the  boundary  value  problem  of  the  form: 

4  *  e  ^  cos  x  +  h2t  -  j  sin  (k  ^  x  +  h  ^  t)J  (a  14 

where  kj  =  and  hj  * ^  «  The  above  expression  satisfies  Eqs.  A 10,  All,  and 

A12  to  good  approximations  if  the  following  conditions  hold: 


«  1 


«  1  h,  « 


(A  15 


*i  fh 

j  "  Tf 


(A16 


Combining  Eqa.  A9  hnd  A 13,  these  results 

1  8 


1-  2  w  yi-  cos(kjX  +  hjt) 


and  upon  substituting  from  Eq.  14: 


1-t-j —  cos(kjX  +  h^t) 


p  ■■  -  \  '  |  “  J  -  J  j 

Cg  ’J'gT  - -  sin[k2x  +  h2t  -  j  sin(kj  x  +  hjt)], 

^  1  -  2  X  COS  (kj  X  +  hj  t) 


[l  +*  -p-  cos  (kjx+hjt)]  sin  jk2x+h2t- j  sin(kjX  +  hjt)J 

2  1  1  ,• 
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tz'dfr  ^8in  >k 


«*  ^ 

| sin  ^<2  x  +  h2 1  -  j  sin(kjX  +  hjt^  +-J - -  sin^k^  +  k^)x  +  (h^  +  h ^)t 


-j  sin(  kj  x  +  h  )j  + -^-i-  sin  jjjk^  -k  ^)x  +  (h2  -  h  j)t  -  j  sin  (k  j  x  +h  jt)J  ^ 


(A  17 


Note  that  these  terms  are  similar  to  the  frequency  modulated  waves  of  radio  systems, 

al 

and  that  the  added  terms  are  linear  in  -r —  .  The  deviation  from  the  first  order  theory 

1 


is  thus  more  serious  when  several  components  are  present  in  the  wave  structure,  since 
the  perturbations  of  a  single  component  due  to  finite  amplitude  are  only-  proportional 


to 


This  large  "  intermodulation  distortion" 


may  not  be  entirely  disadvantageous 


as  far  as  the  detection  of  an  submersed  object  is  concerned,  since  the  moving  object 
may  cause  the  long  waves  which  modify  the  natural  waves  of  much  shorter  wavelength. 


